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Abstract 

Studies of D and B mesons decays into hadrons have been used to test the standard model 
in the last fifteen years. A heavy meson decay involves the combined effects of a primary weak 
vertex and subsequent hadronic final state interactions, which determine the shapes of Dalitz 
plots. The fact that final products involve light mesons indicates that the QCD vacuum is an 
active part of the problem. This makes the description of these processes rather involved and, 
in spite of its importance, phenomenological analyses tend to rely on crude models. Our group 
produced, some time ago, a schematic calculation of the decay D+ —)■ A“7r'''7r''', which provided 
a reasonable description of data. Its main assumption was the dominance of the weak vector- 
current, which yields a non-factorizable interaction. Here we refine that calculation by including 
the correct momentum dependence of the weak vertex and extending the energy ranges of Trvr and 
Ktt subamplitudes present into the problem. These new features make the present treatment more 
realistic and bring theory closer to data. 
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I. MOTIVATION 


Non-perturbative QCD calculations are difficult and can only be performed in approxi¬ 
mate frameworks. The grouping of quarks into two sets, according to their masses, provides 
a convenient point of departure for approximations. Quarks u, d, and s can be considered 
as light and quarks c, b, and t, as heavy, even though the s-quark is not too light and the 
c-quark is not too heavy. This approach is useful because light quark condensates are active 
close to the ground state of QCD and give rise to highly collective interactions. 

Pions and kaons are the most prominent light quark systems, but data available for elastic 
Kn scattering are scarce and decades old. They were obtained from the LASS spectrometer 
at SLAC[1, 2], in the range 0.825 < y/s < 1.960 GeV, by isolating one-pion exchanges in the 
reaction KN — )■ ttKN. In the last ten years, information about Ktt interactions was also 
produced by hadronic decays of D mesons. In particular, data from the E791 and FOCUS 
collaborations[3, 4] for the reaction D+ —)■ allowed the S'-wave Kn sub-amplitude 

to be extracted continuously from threshold up to the high energy border of the Dalitz plot. 
Hope was then raised that these data could improve the description of elastic Kn scattering. 
However, decay data differ signihcantly from those given by the LASS experiment and this 
discrepancy motivates our interest in this problem. 

The description of the decay —)■ K~n~^n~^ must include both the weak c —)■ s vertex 

and hadronic hnal state interactions (FSIs), which correspond to strong processes occur¬ 
ring between primary decay and detection. The study of weak vertices departs from the 
topological structures given by Chao[5], which implement CKM quark mixing for processes 
involving a single W. As primary decays occur in the presence of light quark condensates, 
the direct incorporation of Chao’s scheme into calculations is not trivial and one is forced 
into hadronic descriptions. These include both the use of form factors in weak vertices, 
as in the work of Bauer, Stich and Wirbel[6], and the treatment of relativistic hnal state 
interactions. High-energy few-body calculations begin to be available now[7-9] and several 
works have already employed held theory to FSIs in heavy meson decays[10-17]. 

In this work, the decay —)■ K~n^n^ is treated by means of chiral ehective lagrangians, 
supplemented by phenomenological form factors. This framework is motivated by the small¬ 
ness of the u, d, and s masses, when compared with the QCD scale A ~ 1 GeV. The light 
sector of the theory is therefore not far from the massless limit, which is symmetric under the 
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chiral SU{3) x SU{3) flavour group. In this approach, light condensates arise naturally and 
pseudoscalar mesons are described as Goldstone bosons. Quark masses are incorporated 
perturbatively into effective lagrangians[18, 19], whereas weak interactions are treated as 
external sources. Chiral perturbation theory was originally designed to describe low-energy 
interactions, where it yields the most reliable representation of QCD available at present. 
Its scope was later enlarged, with the inclusion of resonances as chiral corrections[20], and 
the unitary ressummation of diagrams[21]. Suitable coupling schemes also allow the incor¬ 
poration of heavy mesons[22]. 

A similar theoretical framework has already been employed by our group [16], in an ex¬ 
ploratory study of FSIs in —)■ K~7r^7r^. With the purpose of taming an involved 

calculation, in that work we made a number of simplifying assumptions. Among them, the 
weak vertices were taken to be constants, isospin 3/2 and P waves were not included in 
intermediate Ktt amplitudes, and couplings to either vector mesons or inelastic channels 
were neglected. In spite of these limitations, that work allowed the identification of leading 
dynamical mechanisms and gave rise to results which are reasonable for the modulus and 
good for the phase of the S'-wave Kn sub-amplitude[3, 4]. In this work, we focus on the 
vector weak amplitude and improve the description of the weak vertex, by including both 
the correct momentum dependence and better phenomenology for an intermediate tttt sub- 
amplitude, and the description of a Kn subamplitude at higher energies. These new features 
tend to reduce the gap between theory and experiment. 

II. DYNAMICS 

We denote by [K~n~^]s, the A-wave K~n~^ sub-amplitude in the decay D~^ —)■ K~n~^n~^, 
which has been extracted by the E791[3] and FOCUS [4] collaborations. The decay begins 
with the primary quark transition c —)■ s W~^, which is subsequently dressed into hadrons, 
owing to the surrounding light quark condensate. In the absence of form factors, this 
structure gives rise to the colour allowed process shown in Fig.l, where a and b involve an 
axial current and c contains a vector current. As one of the pions in diagram c is neutral, it 
does not contribute at tree level. 

Inclusion of final state interactions, due to successive elastic Kn scatterings[16], yields 
three families of diagrams, as in Figs.2. It is worth noting that these series do not represent a 


3 




FIG. 1: Topologies for the weak vertex: the dotted line is a scalar resonance and the wavy line is 
the W~^, which is contracted to a point in calculations. 

axial 
series 



vector 

series 




FIG. 2: Final state interactions starting from the axial weak vertex (axial series ) and from the 
vector weak vertex (vector series); in the former, the pion plugged to the W~^ is always positive, 
whereas the K inside the loop can be either positive or neutral; in the latter, the tree diagram does 
not contribute, since one of the pions plugged to the W~^ is neutral. 

loop expansion, because loops are also present within the Ktt amplitude. The W~^ is shown 
explicitly, just to indicate the various topologies, and becomes point-like in calculations. 
A family of FSIs endows the forward propagating resonance in Fig. lb with a dynamical 
width[23]. Processes involving resonances have already been considered in Refs.[13, 14, 17], 
whereas quasi two-body axial FSIs were discussed Ref. [15]. An important lesson drawn from 
our previous study[16] is that, for some yet unknown reason, the vector weak amplitude, 
represented by diagram (c) of Fig. 1, seems to be favoured by data[4]. This amplitude receives 
no contribution at tree level, since the emitted by the c-quark decays into a pair. 
Therefore, leading terms in this process necessarily involve loops, which bring imaginary 
components into the amplitude. 

The hrst non-vanishing contribution to the vector series is given in Fig.3. As the W is 
very heavy, one keeps just hadronic propagators, which render loop integrals finite. Denoting 
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FIG. 3: Leading vector current contribution, dressed by form factors and Trvr interactions (in the 
small green blob). 

by Aq the amplitude for the process K^7i^7i~^ without FSIs and by Txn that for 

—)■ 7r~^K~, the amplitude Ai of Fig.3 can be schematically written as 

r 

M = j j^.TLK^kAo, ( 1 ) 

where P is the loop variable and A,r and /Ak are pion and kaon propagators. 

The amplitude Aq is described in App.B. The D ^ WK vertex includes D* intermediate 
states, associated with form factors parametrized in terms of nearest pole dominance [24] 
and could be a vector or a scalar. The W titt form factor is shown in Fig.4 and includes 
the p, with a dynamical width. The bare resonance is treated employing the formalism 
developed in Ref. [20] and its width is constructed using the P-wave elastic tttt amplitude. 
The IF —)■ TTTT form factor is time-like and its inclusion into the vector series of Fig.2 can, 
in principle, give rise to final state interactions depending on both tttt and Ktt amplitudes. 
With the purpose of keeping complications to a minimum, we consider just tttt interactions 
which are contiguous to the W and occur to the left of the hrst Kti amplitude. 



FIG. 4: Structure of the W —)• ttyt form factor; the blue blob is the elastic tt-tt amplitude. 

The evaluation of Fig.3 requires the iFvr amplitude in the interval 0.401 GeV^ < s < 
2.993 GeV^. As LASS data[l] begins only at s = 0.681 GeV^, one covers the low-energy 
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region by means of theoretical amplitudes, based on unitarized chiral symmetry[20]. Our 
intermediate S'-wave Kn amplitude, denoted by is thoroughly discussed in App.C. 
Using results {B17) into eq.(l), one hnds 



[Gf cos^ecF^'^iO) 



d'^i 1 ml 

(27r)4 D^Dk Dp 
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2iVf^ + — 2 ?ti^2 
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( 2 ) 


where Gp is the Fermi constant, Oq is the Cabibbo angle, is a coupling constant[24], 

the factor \/2/3 is associated with the transition K~n^, whereas = [(^—pa)^ — 

M^j, Dk = [{i — Py — M]^], Dy = [i’^ — my], Ds = W' — mW, in which the subscripts V and 
S stand for the D*(2112) and Zi)*Q(2317) states. Finally, Dp is a complex function dehned 
by eqs.(B15) and (B16). This structure yields 


= -i a 
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IOtt^ 


{ P ^wKpV ~ ^nKV + 


wpv F ^KpV 7 IrpVs} ) ( 3 ) 


with 


a = [Gp cos^ 6*c F’/^'^(0)] mlmy , 
/3 = M‘l) + 2Ml + M‘^-ml-2ml2 , 
7 = - {Ml-M]^) {ml-ml)/ml , 


and 


jS — 
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16n^ 
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-‘■abed 


(27r)4 DaDbDc ’ 

The form of these integrals is discussed in App.D. 


d^i IOtt^ 

(27r)^ DaDbDcDd ' 


(4) 

(5) 

( 6 ) 

(7) 


III. VECTOR ESI SERIES 

In the decay D-^ —)■ iF“7r’''7r’'', there is no tree contribution to the vector FSI series, as 
in Fig.2. However, before moving into this reaction, it is instructive to assess the relative 
importance of allowed tree and one-loop contributions in the decay D"*" —)■ indicated 
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FIG. 5: Vector current diagrams contributing to the decay 




TT vr' 


in Fig.5. The amplitude describing the left diagram is denoted by Aq and given in eq.(B17). 
Projecting out the S'-wave, we hnd 


A 
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2(3 [ml; e]-al^+(3 


( 8 ) 

(9) 


where 6i are complex parameters given in table I (App.B). The hrst order amplitude is 
obtained by replacing the isospin factor \/2/3 with —1/3 in eq.(3) and reads 


= so 








s 

TvKpV ^ttKV 




+ hpV + ^ 


KpV 


-^t^KpVs] 


( 10 ) 


Results for the moduli of Ag and Af, displayed in Fig.6, indicate a clear dominance of the 
former. The main structural difference between both terms is the factor |T^^/487r^} in the 
latter, associated with a hnal state scattering. Its scale can be understood by noting that 
chiral symmetry predicts this amplitude to be = 2 M^^Mk/F'^ ~ 13 at threshold whereas 

LASS data [1] indicate that it reaches a maximum of ~ 50 around mi 2 ~ 1.33 GeV. 
Therefore the factor |T^^/487r^} is always smaller than 1/10 and pushes down the loop 
contribution. This result can be taken as an indication that the vector series, as given in 
£g.2, converges rapidly. The conhrmation of this hint depends, of course, on the explicit 
calculation of next terms in the series. 
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FIG. 6: Modulus of the D~^ —>• amplitude (full line) and partial contributions from eqs.(8) 

(dashed line) and (10) (dotted line). 


IV. RESULTS - S-WAYB 

One of the purposes of this work is to understand the role played by the high energy 
components of intermediate tttt and Ktt subamplitudes in the description of data. Predic¬ 
tions from eq.(3) for the phase and modulus of [K~7i~^]s, the S'-wave K~7i~^ sub-amplitude 
in D~^ —)■ are given in Figs.7 and 8. 

As far as the tttt subsystem is concerned, the data of Hyams et ah [25] are used in a 
parametrized form, in the whole region of interest, as discussed in App.B. For the sake of 
producing a contrast, we also show curves corresponding to the low-energy vector-meson- 
dominance approximation, in which the P-wave amplitude is described by just an interme¬ 
diate p-meson, which amounts to using just the hrst term in eq.(B16). In the case of the 
Kit amplitude, data are not available for energies below 0.825 GeV [1] and two alternative 
extensions are given in App.C. One of them is based on a two-resonance £t, which encom¬ 
passes both low- and high-energy sectors, whereas in the other one LASS data[l] is used 
directly, when available, and extrapolated to the threshold region by means of a £t. In the 
sequence we refer to these versions as fitted and hybrid, respectively. The main difference 
between them is that the former excludes points around P ~ 1.7 GeV, shown in Fig. 11, 
where two-body unitarity is violated. 

Inspecting the hgures, one learns that the improvement in tttt phenomenology is more im- 
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FIG. 7: Predictions for ^ K ~phase (full blue curve), based on the parametrized Trvr and 
Ktt amplitudes given in appendices B and C, compared with FOCUS data[4]; the blue dotted curve 
is the previous one shifted by —55*^; the dashed blue curve is based on one-p pole approximation 
for the vrvr amplitude; in the red symbol-continuous curve the hybrid model was used for the Kn 
amplitude. 

portant for the modulus, where it influences considerably the curve behaviour and increases 
signihcantly the range in energy where the theoretical description proves to be reasonable. 
In the case of the phase, effects associated with tttt phenomenology are small and visible 
only above mi 2 > 1.2 GeV. On the other hand, the use of either the htted or hybrid Kn 
amplitudes produces equivalent results, except at the high energy end, where none of them is 
satisfactory. This seems to indicate missing structures, that could be associated with other 
topologies in ^ K~7r~^7r~^ decay. 

As experimental results for the FOCUS phase[4] include an arbitrary constant, in Fig.7 
we also show our main result displaced by —55*^. One notices an overall good agreement with 
data, from threshold up to mi 2 ~ 1.4 GeV. As our results were based on the vector series 
shown in Fig.2, which does not contain a tree contribution, there are two sources of complex 
phases in this problem. One of them is that associated with the Kn amplitude, whereas the 
other one is less usual and due to the loop including the weak vertex. Our results indicate 
that the latter is rather important over the whole energy range considered. This shows the 
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FIG. 8: Predictions for D'^ ^ K ~modulus (full blue curve), based on the parametrized ttyt 
and Ktt amplitudes given in appendices B and C, compared with FOCUS data[4], using arbitrary 
normalization; the dashed blue curve is based on one-p pole approximation for the vrvr amplitude; 
in the red symbol-continuous curve the hybrid model was used for the Ktt amplitude. 

relevance of proper three body interactions, which share the initial momentnm with all hnal 
particles at once. 


V. CONCLUSIONS 

In this work we calcnlate the weak vector cnrrent contribntion to the process —>■ 

, employing intermediate tttt and Ktt intermediate snb-amplitndes valid within 
most of the Dalitz plot. Together with the nse of a proper P—wave weak vertex, this 
extends a previons stndy made on the snbject[16]. We still concentrate on [K~tt^]s.i the 
S'-wave K~tt^ snb-amplitnde, and present predictions for both the phase and modnlns, 
given by the bine cnrves in Figs.7 and 8, are qnite satisfactory from threshold to 1.4 GeV. 
Resnlts for the modnlns, in particular, improve considerably our previous hndings, showing 
that intermediate tttt subamplitudes are important and need to be treated carefully. As 
far as the phase is concerned, the most prominent feature is the fact that it has a large 
negative value at threshold. In QCD, loops are the only source of complex amplitudes and, 
in this problem, the energy available in the loop of Fig.3 can be larger than both Ktt and 
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Kp thresholds. This yields a rich complex structure for the loop containing the W, with a 
phase ©L which adds to the phase Qk-k already present in the intermediate Kn amplitude. 
Therefore Ql represents the gap between the two and three-body phases, which depends on 
both mi 2 and 77123 , showing that Watson’s theorem does not apply to this case. Our results 
both conhrm the dominance of weak vector currents in this branch of decays and indicate 
that proper three body hnal-state interactions, in which the initial four-momentum of the 
D~^ is shared among all hnal particles, are rather important over the whole energy range 
considered. In a parallel study, to be presented elsewhere, we found that this feature is also 
present in the P—wave projection of hnal-state Kn subamplitude, which has a non-vanishing 
phase at threshold. 
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Appendix A: kinematics 


The momentum of the P-meson is P, whereas those of the outgoing kaon and pions are 
P 2 , Pi, and ps , respectively. The invariant masses read 


^12 = {P 1 +P 2 Y = Ml + -h 2pi-p2 , 

(Al) 

= (Pi +Pzf = 2 + 2 Pi -ps , 

(A2) 

ml'i = {P2+Plif = Ml + Ml + 2p2-p3 , 

(A3) 

and satisfy the constraint 


Ml = + ^13 + ml^-2Ml- Ml . 

(A4) 

The projection into partial waves for subsystem (12) is performed by going to its center of 

mass and writing 

^13 = «i 3 - f^i 2 cos 9 , 

(AS) 

^23 — ^23 + /^12 COS 9 , 

(A6) 
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053= [Ml + 2Ml + Ml-ml,- (Mj-Mj) (Mj-Mjj/mfJ/2 , (A7) 

a'23 = IK + 2Ml + Ml - ml, + (Ml-Ml) (Ml-Ml)lml,] /2 , (A8) 

ft2 = 2l5V. (A9) 

q' = [,< - 2(Ml+Ml) ml, + (Ml-Mlr^] , (AlO) 

/ Y 771^2 

Q' = ^ [m^2 -2(M2 + M|,)m?2 + , (AH) 

2A/m|2 

where 0 is the angle between the momenta of the pions. 


Appendix B: basic D+ —)■ A°7r°7r+ amplitude 


Our description of the decay —)■ A'“7r’''7r’'' includes both the primary weak vertex 

and hadronic hnal state interactions, associated with successive Ktt scatterings. When the 
W —)■ TTTT vertex is corrected by means of time-like form factors, both the p-meson and 
P-wave TTTT interactions also become part of the problem. This could, in principle, give 
rise to a structure of hnal interactions depending on both tttt and Ktt amplitudes. For the 
sake of keeping complications under control, we consider here just tttt interactions which 
occur to the left of Ktt amplitudes. Therefore, the amplitude for the process D~^(P) —)■ 
-A°(pA')7r°(po)7J'^(p+), given in Fig.4 and denoted by Aq, becomes the basic building block 
in the evaluation of the weak vector series. 

We begin by constructing the isospin / = 1, P-wave tttt amplitude. The momenta 
of the outgoing pions are p+ and po , whereas those inside the two-pion loop are g+ and 
go- The total momentum is Q = (p+ + Po) = (q+ + Qo) and the loop integration variable is 
£ = (g+ — q -)/ 2 . Assuming that, at low energies, tttt interactions are dominated by a 0 { q ^) 
contact term supplemented by the 0 { q '^) p-pole contribution, the effective lagrangians in 
Ref. [20] yield the tree contribution 


= it 


u] 


1 


2Gl 


V 


F‘^ s — rn?p 


(Bl) 


where P is the pion decay constant and Gy describes the p tttt coupling. The approximation 
Gy = F/V2 66 MeV yields a more compact structure, given by 


pi = 


P2 


m: 


s — mi 


(B2) 
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For free particles in the center of mass frame, {t—u) = (s—4M^) cos 6 and P-wave projection 
yields the kernel 


j^pi ^ _ (f—^33^ 

3 P^ s — m?p 

The iteration of this kernel by means of intermediate two-pion states produces the uni- 
tarized amplitude [21] 


rjiPl _ 

TTTT 




where is a divergent loop function. Therefore, we write it as the sum of an inhnite 
constant Aoo and a regular component 0 ,^ 71 , given by[16] 

A s Ta , f!,-2Ay + VAl , . Va] 

- 2Ml - ’ 


\ = s^-AsMl , 


where P = 1/2 is the symmetry factor for identical particles. After regularization, one hnds 


rjiPl _ 

TTTT 


1 + (^TrTr + IPtitt) 


where is an arbitrary constant. This amplitude is related with phase shifts by 

r 1 

(B7) 

and we fix by the phase at 90°. The 1=1 amplitude to be used in the evaluation of 
Aq is given by eq.(B6) multiplied by (3 cos6'). It is denoted by and can be cast in the 
covariant form 

rj .1 _ o - u) rT.Pl _ ^ iP+ - PoT^U rT.Pl ARoA 

s-AMl S-4M2 ^ 

Going back to the decay amplitude Aq and reading the diagrams in Fig.4, one hnds 


r 2 

— Tfl 

ATT = [GFCos'^ec]{K^\V^\D^) — -^4 

P ?Tl„ 


1 .A ,PI 


(p+-Po)m + *6 


— 4M^ W ) 1 


i H.U - 


{2nY [(£+g/2)2 - Ml] [{^-Q/2f - Ml] ’ 
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where Gp is the Fermi constant, 6c is the Cabibbo angle, is the weak vector current. 
The regular part of can be related with eq.(B5) and one has[26, 27] 


v = + (BIO) 

Using this result into eq.(B9) and recalling that (p+ —poYQu = 0 for on shell particles, one 
has 

9 

— 7Tl 

Ao = [Gf cos^^c] ( K'^IV^ID^ ) {p+-po), -=r ’ 

Dp = {Q^-mD - (mJ/3F2) {Q^-4M^) + a.] • (B12) 


The vector current matrix element is written as 


{K^\V^\D-^ ) = iP^+ [FF^{Q") - F.^^iQ")] (B13) 


and form factors are parametrized in terms of vector and scalar cs nearest poles as [24] 


FiiQ^) = 


1 - 


and id)(g^) = 


1 — Q'^/rrFg 


(B14) 


with V D;(2100), S T);(2317) and F^^(O) = 0.75. 

The denominator Dp describes the p meson and includes its dynamically generated width. 
The function Dp does not vanish along the real axis, in spite of the bare p propagators in 
Fig.4. It has a zero in the second Riemann sheet, quite close to the value quoted in Ref.[28], 
namely at = {rrip — i Fp/2)^ , rrip = 762.4 ± 1.8 MeV, Fp = 145.2 ± 2.8 MeV. 

In order to simplify calculations one notes that the ratio w?p/Dp in eq.(Bll) is related to 
the P-wave amplitude given by eq.(B4) by 


" V _ _ ^ rjnPl 

Dp~ s- 4M2 ■ 

Using the data from Hyams et al.[25], we htted this ratio using the structure 

ml _ Np N, N2 

Dp s — 6p s — 9i s — 6*2 ’ 


(BIS) 


(B16) 


with free parameters 6^ = dkR + i 9ki and + i N^j , given in table I. 

In Figs. 9 and 10 we display the importance of the inclusion of higher poles in eq.(B16) 
in extending the agreement with Hyams et. ah[25] data. 
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k 


0i 

Nr 

Nj 

P 

0.580133 

-0.1137172 

0.6131598 

-0.1107509 

1 

2.539625 

-0.6468928 

0.2418401 

-0.1080483 

2 

3.642091 

-0.1595399 

0.0016668 

-0.1941643 


TABLE I: Fitted parameters in eq.(B16) 




FIG. 9: Results for vrvr phase and modulo with only one-p (dashed) and adding another two poles 
(continuous), compared with Hyams et al.[25] (dotted). 

The expression for Aq to be used in calculations is obtained by assembling previous 
results, and one finds 
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FIG. 10: Function m^p/Dp] the red continuous curve represents eq.(B15) with parameters from 
Hyams et al.[25] and the black dotted curve is our fit, using eq.(B16); as data begin at 0.6GeV, 
the red curve to the left of the vertical dashed line corresponds to an extrapolation. 

Appendix C: Ktt amplitude 


In this work, one needs the elastic Kn amplitude over the full Dalitz plot. As there are 
no data[l] available in the interval 0.401 GeV^ < s < 0.681 GeV^, one encompasses this 
region with the help of a theoretical amplitude, based on the unitarized chiral symmetry. 
This model has been discussed in detail in Ref. [16, 29] and here we just summarize its main 
features. 

For each spin-isospin channel, the unitary amplitude Tli is obtained by ressumming the 
infinite geometric series 


Tli = K-li — ICli JCli + K-li [fiKn+CLi] JCli [^k-k+Cli] K-li + ■ ■ ■ 

K-li 


(Cl) 


1 + [^Ktt+Cli] ICli 

where Kli is a kernel and the function CLk-ki related with the two-meson propagator, is given 
by[16] 


M^-Ml Mk 


M. 


K 


x/A 


In 






+ 171 - 
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A = s' - 2s(M^ + M|) + {Ml-Mlf , 


(C2) 


and Cu is a constant. Chiral perturbation theory determines the kernels ICn as the sum 
of a 0{q^) contact term[19], supplemented by 0{q^) corrections, which we assume to be 
dominated by s-, t- and w-channel resonances[20]. In order to £t LASS data[l], we also 
included a higher mass resonance, as described in Ref. [29]. 

In the case of the S' 1/2 wave (L, / = 0 , 1 / 2 ), the theoretical kernel is written as = 

fCsG + , where /Cbg is a real background and fCu includes resonances. The former is 

given by /Cbg = + /C 5 + cy /Cy , with 

9 V 2 - (Ml+Ml)] , (C3) 

^s = - 2cd {cd-Cm){Ml + Ml) - 2c], q^] 

+ [cdml-2{cd-Cra)Ml] [cdml-2{cd-Cm)M‘j^] 4(q';mo)} 

+ ^ - ‘2cd{cd-c^){Ml + Ml) - 2c]q^] 

+ [cdwl - 2{cd-Cra)Ml] [qTT ig - 2 (q-c,„)M|] 

+ { '^d(^^*+2M/ + M|;-s + 2q') + 2crf(Q-Cm) 

+ Cdm\.-{cd-cJ) {Ml + Ml) , (C4) 



(C5) 

(C6) 

(C7) 
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where F, q, Cm, q, Cm, and Gy are coupling constants and the CM three-momentum is 


= — [s" - 2s (A^+Mj) + (Ml-Mlf] 


(C8) 


Two s-channel resonances are incorporated as sum of Breit-Wigner functions[29] 


Kh = 


2F4 


[CdS- {Cd-Cm) {MI + M]^)Y [CdhS- jCdb-Cmb) + 
s-mj^.+iglQa/8^^^/s s - ml + i gf Qb/STiy/s 


9i -^i T ) 

Qi 


^ (1 - ht/s) . 


, (C9) 

(CIO) 

(Cll) 


The usual inelasticity parameter p, evaluated for S' 1/2 data, is shown in Fig.ll. Points for 
which p > 1 within error bars were discarded in our £t. 



FIG. 11; Inelasticity parameter g for S 1/2 LASS data. 

We have extended S' 1/2 Ktt data to threshold by means of two different fits. The first one 
includes just a single resonance and holds for energies smaller than 1.33GeV, whereas the 
second one includes two resonances and is valid over the whole Dalitz plot. Their correspond 
respectively to x^/n.d./. = 0.55 and y^/n.d./. = 1.62. Our parameters, in suitable powers 
of GeV, are: F = 1.02722, Gy = 0.0686287,m8 = mo = 0.983 and = 1.124899 x lO-^, 
mx* = 1.108858, q = 0.0254505, Cm = 0.1483455, Aa = 4.563646, Ba = -2.055842, ha = 
1.138489, cy = 0.26200 for the single resonance fit and G^^^^ = —2.273182 x 10“^, mx* = 
1.338404, Cd = 0.026607, = 0.017428, Aa = 4.952313, Ba = -1.956429, ha = 1.130126, 

mb = 2.003338, Cdb = 0, Cmb = 0.166268, Ab = 5.042537, Bb = —7.182061, hb = 1.809129, and 
Cy = 0.89272, for the two-resonance case. 


18 




modulus 


Both fits for the modulus an phase are given in Fig.12. In the —)■ decay 

amplitude, alternatively, we can use directly empirical data from LASS[1] and merge it with 
the low energy £t, where there is no data. This became what we called hybrid amplitude. 




FIG. 12: Fits for the modulus and phase of the K-k S' 1/2 LASS data; points within the regions 
indicated as gap in the top axis were excluded from the fit. 

In Fig. 13 we show the real and imaginary components of the amplitude. One notice 
that values for the real part at threshold are different, namely 24 and 30, and they can be 
compared with those obtained by ChPT[31] and dispersion relations[30], respectively T = 
21.7 and T = 25.5. These values indicate that the low energy £t is more suitable to describe 
low energy behaviour. It is worth noting that both real and imaginary components are very 



FIG. 13: Real and imaginary components of the S 1/2 Ktt amplitude fitted to LASS data (squares) 
and extended to low-energies using chiral symmetry. 

small for ^/s ~ 1.65 GeV in the two-resonance result. 
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Appendix D: loop integrals 


/■S' _ 

abc 


dH 




(27r)4 DaDbDc 


[S — 

abed 


dH 




(27r)4 DaDbD^Dd 


We begin by discussing the integrals , given by eqs.(7). Their treatment can be sim- 
plihed because the p and the cs state entering the form factor share the same momentum. 
This allows, for instance, one to write 


InKpV [4x1/ 4xJ , (Dl) 

where 0 is the parameter dehned in App.B. Similar simplihcations can be performed every 
time subscripts pV or pS occur. 

The integral is important in this problem because its imaginary part is determined 
by two different thresholds, associated with cuts along Kn and Kp propagators. Using 
results from App.B, one writes 


/■S' ^ 

^ 7 tKp 


dH 




Nn 


(D2) 


(27r)4 [(£-p3)2_M2][(£-P)2-Mi] + ■ 

Representing this function by means of Feynman parameters and performing one of the 
integrals analytically, one hnds 


— iN^U 


p ^^irKp ? 


with 


and 


hl-TrXp / det d,rXp(^) 5 

Jo 


J-KKpifl^ 

Oj 

Fip 


1 

71 


In 


141 

IG 2 I 

Gi 

1141 


+ 44 


Gig = 


— [ Opi — 0 f 2 — dci + 0g2 ] 

|2Mj,a + aT Va| i,„ , 
Ml 

\B q: \/A| 


Ml 


(D3) 

(D4) 

(D5) 

(D6) 

(D7) 

(D8) 


B = [Sr - sS; - - Ml + - a (Mj - Ml + mjj)] , (D9) 
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X = -AMlC , (DIO) 

C = [(1 —a) + a — a (1 —a) m ^2 ] • (DU) 


The p width is incorporated into the factors Np, 0/, and the case of a point-like resonance 
is recovered by making iVp —)■ 1,0/ —)■ e. 

The integral I^xv 

d^^ IOtt^ 1 


r* — 

^ttKV — 


= m 


(27r)4 [(£-p 3)2-M2] [(£-P)2_Mi] £2 _ ^2^ ’ (^^2) 

and its evaluation is totally similar. However, as now my > , its imaginary part comes 

just from the cut of the diagram along the Kn subsystem. Integrals I^py, if^pv and 
not depend on 771^3 • 
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